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STRONG  INTERACTION  ON  A  PLATE  WITH  CONSIDERATION 
OF  SLIP  AND  NEAR-WALL  TEMPERATURE  JUMP 

A.  A.  Bogacheva,  V.  S.  Galkin 

We  will  solve  the  problem  of  the  motion  of  a  perfect  gas  In  a 
laminar  boundary  layer  on  a  semi-infinite  flat  plate  at  zero  angle 
of  attack,  taking  Into  account  the  slip  and  jump  of  gas  temperature 
at  the  wall  and  the  strong  interaction  between  the  boundary  layer 
and  a  hypersonic  lnviscid  flow.  It  is  assumed  that  the  viscosity 
coefficient  u  linearly  depends  on  temperature  T,  Prandtl  number  Pr  ■  1. 
In  order  to  calculate  pressure  P  we  used  the  method  of  tangent  wedges. 
The  slip  and  near-wall  temperature  Jump  is  taken  into  account  by 
linearization  of  the  solution  of  the  relatively  well-known  self-similar 
solution  with  the  usual  boundary  conditions  of  "adhesion"  *  (no  slip)  . 
We  will  examine  the  case  of  an  insulated  plate  and  the  case  of  a  cooled 
(Tw  <  T0)  plate  (plate  temperature  Tw  =  constant.  To  is  stagnation 
temperature  of  the  free-stream  flow,  i.e.,  the  temperature  of  the 
insulated  plate) . 

*  We  are  concerned  with  a  self-similar  solution  for  a  strong 
interaction  of  zero  order  in  the  nomenclature  of  Hayes*  work  [7]. 
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This  statement  of  the  problem  was  used  earlier  in  Shidlovskly's 
work  [1]  for  the  case  Tw  =  T0.  Calculations  were  conducted  at  a  value 
of  the  specific  heat  ratio  x  *  1.4  and  x  -  5/3,  and  the  solution  of 
the  zero-approximation  equation  was  found  by  interpolation  of  the  well- 
known  Falkner-Skan  solution  tables  for  other  values  of  the  parameter 
entering  Into  this  equation.  Obviously,  mistakes  were  made  during 
this  interpolation  which  led  to  qualitatively  false  conclusions:  slip 
Increases  pressure  P  and  the  boundary-layer  displacement  thickness  6, 
whereas  the  nature  of  its  effect  on  the  local  coefficient  of  frictional 
drag  Cf  depends  on  the  magnitude  of  x.  In  this  work  equations  of  the 
zero  and  first  approximations  were  solved  simultaneously  in  the  elec¬ 
tronic  computer  which,  incidentally,  allowed  us  to  refine  the  data  of 
works  [2],  [5]  for  a  strong  interaction  without  regard  to  the  rare¬ 
faction  factors.  The  results  of  the  calculations  for  x  *  1.4  and 
x  =  5/5  show  that  when  Tw  £  To*  slip  and  temperature  Junq>  decrease  P, 

6,  and  especially  the  heat-transfer  coefficient  from  the  gas  to  the 
plate . 

1.  We  will  examine  motion  in  a  laminar  boundary  layer  on  a  plate 
under  conditions  of  strong  interaction  of  the  boundary  layer  with  a 
inviscid  flow.  A  strong  interaction  occurs  [4]  when  (dfi/dx)a  >  1, 
where  is  the  free-3tream  Mach  number,  x  is  the  coordinate  along  the 
plate.  In  addition,  magnitudes  of  the  order  must  be  negligibly 
small  in  comparison  with  unity  (we  assume  that  the  magnitude  2/(x  —  1) 
is  not  one  of  the  determining  parameters  when  evaluating  orders  of 
magnitudes,  i.e.,  2/(x  —  1)  ~  l)  .  Otherwise  strong  interaction 
occurs  at  the  leading  edge  of  the  plate  being  passed  by  a  flow  with 
»  1,  where  the  rarefaction  effects,  caused  by  slip  and  gas  tem¬ 
perature  Jump  on  the  plate,  and  the  entropy  effects,  caused  by 
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vortlcity  of  the  inviscid  flow  on  the  outer  edge  of  the  boundary- 
layer  flow,  can  be  particularly  significant. 

The  Investigation  of  these  effects  can  be  facilitated  since  we 
can  study  them  individually  by  means  of  boundary-layer  equations,  [1], 
[5]-[7]j  while  appropriately  changing  the  boundary  conditions,  since 
the  effects  of  rarefaction  and  vorticlty  are  small,  approximately  of 
order  one.  Actually,  let  the  terms  of  the  boudary-layer  equation  be 
of  the  order  of  unity.  Then  these  rarefaction  effects  are  maximum 
in  a  case  of  strongly  heated  plate  [5]  when  they  are  of  the  order  of 
6/x.  The  order  of  magnitude  of  entropy  effects  is  equal  [7]  to 

(6/x)*-V(**) . 


In  the  zero  approximation  magnitudes  of  the  order  of  6/x  and 
higher  are  disregarded.  Here,  the  boundary-layer  equations  have  a 
self -similar  solution.  Terms  of  the  order  of  6/x  are  taken  into 
account  in  the  first  approximation.  By  virtue  of  their  smallness  we 
can  linearize  the  problem  with  regard  to  the  zero  approximation, 
i.e.,  we  can  investigate  individually  the  effects  of  rarefaction  and 
vortlcity,  while  appropriately  changing  the  boundary  conditions  for 
boundary-layer  equations. 

Therefore,  when  investigating  the  effects  of  rarefaction  the 
boundary  conditions  on  the  outer  edge  of  the  boundary  layer  remain 
the  same  as  in  the  zero  approximation:  the  velocity  at  the  edge  of 
the  boundary  layer  is  u^  *»  uw  total  enthalpy  **  u£/2,  i.e.,  the 
temperature  T^  «*  0.  The  boundary  conditions  of  slip  must  be  fulfilled 
on  the  plate  [8] 


AT  »  T  /  *  -  Tw  =  b  /2 1=2.  ( JL  \ 

**  Y  2  9  \  pa  dy  )f 


/«X  2 —  q  /  ]|  £u\  , 

2  3  \?a  dy)yJ 


2  x  2  —  «  a  * 


(1.1) 


K  +  l 


2-0  ' 
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where  a  Is  the  speed  of  sound,  p  Is  the  density,  a  is  the  accomodation 
coefficient,  a  is  the  reflection  coefficient.  The  term  (3/*)  (u/pT) 
(dT/3x) ,  negligible  within  the  limits  of  applicability  of  the  bound¬ 
ary-layer  equations  is  discarded  in  the  formula  for  slip  velocity. 

When  solving  the  problem  we  will  look  for  the  dependence  of 
pressure  on  x  in  the  form  of  P*  «  P/P^  “  +  dx-1/*),  since 

in  the  zero  approximation  P*  ~  x"1/2,  and  the  correction,  caused  by 
raref action  is  of  the  order  of  6/x  <**  x”1/4,  Here  coefficient  Po  and 
small  parameter  d  are  subject  to  determination,  Xi  “  X  when  x  ■  1, 

X  -  Re"1/2,  Rem  =  u^x/v,  v  =  We  will  introduce  now  (gener¬ 

alizing  [1])  the  dimensionless  dependent  and  independent  variables  by 
the  following  formulas: 


l  =4P„X1jc‘/' 


(1.2) 


We  will  note  that  the  self-similar  motions  depend  only  on  X. 
Using  variables  (1.2)  in  boundary  conditions  (l.l),  when  X  -  0  we 
obtain 


•V*  -  OM'xx. 


AC  =  i>OCx. 


(1.3) 


where,  obviously,  *  is  the  small  parameter  of  the  problem  under  con¬ 
sideration.  Using  further  variables  (1.2)  in  boundary-layer  equation, 
we  will  reduce  them  by  the  usual  method  for  boundary-layer  theory,  to 


jjj  PW'xi  -  n  -  25'F^xx)  =  P  ( 1  +  Yx~,,t)  + 

+  (j-f  j'VuO  +  dr-^ 

2s  [Ci'^x-Cx'Ptl  -  (1+  dx-v.)Cxx- 
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Here  the  terms  of  the  order  of  d2  and  higher  are  discarded. 

In  order  to-  obtain  a  system  of  zero-and  first -approximation 
equations  from  system  (1.4),  we  will  introduce  the  coefficient  k  by 
the  relationship  d-c"1/4  »  ki;  moreover  we  set 

V(i.  0)  -  (~)  ,I‘P(>0  -I-  0;W).  G(>.,  0)  =  St(K)  +  0^,(1). 

Substlutlng  these  formulas  into  relations  (1.3),  (1.4)  discarding 
at  first  terms  of  the  order  of  4,  and  then  of  the  order  of  **  as  com¬ 
pared  to  unity,  we  obtain  a  system  of  equations  and  boundary  con¬ 
ditions  corresponding  to  them  in  zero  and  first  approximations. 

<P"  +  9<P*  +  P(go  —  <p”)  -  0.  gt  f  <t g’t »  0. 

9  (0)  =  <P'  (0)  ■=  0,  gt  (0)  =  TjTt,  <f(oo)=g(oo)  =  1; 

r  +  <p;#  +  (i — 2j5)«p't' + p ft  =  k  [W  + 1  (g'-v*)  ] . 

g\  +  -  —  kg': 

Z  (0)  —  0,  C'(0)  =  9*(0).  gi(0)^bgl(0),  l'(oo)=  gx{oa)  =  0. 

— .  r0/reo-^l/M*e.  -1  =(  )'• 

x  2 

In  order  to  solve  system  (1.6)  it  i3  necessary  to  express  k  in 
terms  of  unknown  functions.  For  this  we  use  formula  P*  =  x  [  (x  +  l)/2] 
M2  (dfi/dx)2,  yielded  by  the  method  of  tangent  wedges  where  the  bound¬ 
ary-layer  displacement  thickness  is 

•  • 

(physically  this  method  of  calculating  pressure  is  Justified  because 
when  calculating  an  inviscid  flow  with  an  error  of  approximately 
(6/x)2  we  can  disregard  flow  through  the  layer,  the  boundary- layer 
thickness  can  be  Identified  with  the  displacement  thickness,  and  we 
can  consider  the  outer  edge  of  the  boundary  layer  as  a  solid  wall. 


(1.5) 

(1.6) 
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Using  these  formulas  we  will  find 


7-*-£{' wT<,'-‘4 

t/5?5'— (1.T) 


A  ™  $(£»  —  <p'*)di.  A  =  5  (ffi  —  tyftdx. 

•  t 


2.  In  order  to  simplify  the  problem  we  will  assume,  generalizing 
[1],  that 

* 

C  *=*  9*  +  *  $  *dX.  ft  ”  ft  +  *jf- 


By  means  of  (1.5)  we  will  reduce  relation  (1.6)  to  the  form 


2*  +  <pz'  +  (1—  23)<p'2  =  — Pg  -f  y(£,— O-Mxp*. 


(2.1) 


S'*  +  (<P«)'  -  —  ft. 


2(0)  =  2(00)  =  0.  g  (oo)  =  0,  £(0)  =  ^=i*;(0). 


(2.2) 


Taking  (1.5)  into  account  we  have  the  following  solution  of  Eq. 

(2.2)* 


g  =  s', (-  X  +  +  ig.(0) +S1  (0)1  gl  { Y  : 


(2.3) 


Here 


*  =  -  4  (TjTt)  (7/.  -  4/.J-,  —  2<p'i)  dX. 

In  order  to  simplify  I2  we  will  integrate  termwise  Eq.  (2.1). 


Then 


p/,  -  <P*  (0)  -r  /  (0)  —  0/2)  /.  -  lim  (92). 


(2.») 
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Prom  the  obvious  finite  condition  of  I2  and  z •  (o)  the  magnitude 
e  =»  lim  (<pz).m\ist  also  be  finM.te.  This  is  especially  clear  in  a  case 
of  an  insulated  plate  when 

«=•©,  /,  -  —  2^'zdX, 

* 

therefore  for  the  existence  oflT  this  integral  it  la  necessary  that 

*»* 

where  m  is  a  certain  constant.. 

In  order  to  determine  thes  value  of  o_  and  the  magnitude  of  g*  (0) 
we  will  investigate  the  conduct  of  functions  £  in  the  vicinity  of 
X  =  and  we  will  examine  onMly  the  principal  terms  of  the  correspond¬ 
ing  asymptotic  expansions.  SJBmce 

So  0)  =  g»  (0)  exp  ^ —  j  <p dX 

then  the  principal  terms  of  thne  items  of  function  £  are 

*i«±£<5.  Q m const.-  (2.5) 

L  2  J  K—m 

Let  us  substitute  exprsssion  (2.3)  instead  of  £  into  the  right- 
hand  side  of  Eq.  (2.1)  .  Then*,  when  9  =  X  -  m,  9'  =  1  this  equation 
with  the  usual  substiutions  far  3uch  an  equation  reduces  to  a  Whittaker 
equation.  The  principal  terac  of  solution  of  a  corresponding  similar 
equation  are 

Dt  (X  -  m)-' Wxp  [-  .  D,(X  — 

the  principal  terms  of  a  partnal  solution  of  a  nonsimilar  equation  are 

where  Dt  and  D2  are  arbitrary  constants,  A  and  n  are  known  constants. 

It  follows  from  these  exapresslons  that  the  solution  is  not  unique 
since  the  boundary  conditions  when  X  =  »  are  satisfied  at  any  values 
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1 


of  g'  (0)  .  It  follows,  however,  that  from  physical  considerations  the 
thickness  of  total  enthalpy  loss  be  finite 

CO  m 

A  =  5  (“/«»)( >  —  C)dX  =  B  —  k(>\g<p'dk, 

*  * 

where  B  Is  the  finite  magnitude.  It  follows  from  this,  with  con¬ 
sideration  of  (2.5)/  that  the  magnitude  of  A  will  be  finite  only  when 
gi  (0)  ■  -g£  (0) .  Assuming  g*  (o)  -  -g£  (0)  and  taking  into  account 
that  Dg  -  0,  follows  from  the  finite  condition  £  we  obtains  e  -  0, 
functions  £  and  £  when  X  -»  »  tend  toward  zero  with  respect  to  the 
exponent  and  Eq.  (2.1)  as  the  final  result  takes  the  form 


z*  +  4*'  +  (1  —  2?)  <p'z  =  (g,  —  9")  +  W’  + 

At  last,  after  several  transformations,  taking  into  account 
formulas  (1.7),  we  obtain 


a /*  a./xACo  -  /x/Mco).  r  =  p,x(\  — p, yi/M «,). 

c,  =  (2/p^D^du/dj/),^  «  C/.X,/,/M-’(I  -C,yx/Mm). 

St  -  ( kdT/dy —  7*)]“*  —  St,  Y% /A4«). 

Q.=  //VP*( 0).  St,  //V4^(0)(l  -  TjTJT1. 

a,  -  - 3/4fcr.  p,  -  -*T.  A  -  - 4.3(9?/. -  4z'(0)  -  4<p'  (0 ^TjT^. 


2  —  0,  /*  ax(x  — 1)  T. 

— * — fT* 


Cf. 


9  &/. 

4  9*(0) 


frr. 


st,  =  2r 


92(0) 

4(4  * 


When  estimating  the  limits  of  applicability  of  these  formulas 
we  must  remember  that  the  small  parameter  of  the  problem  is  ♦  « 
“7^x7^-  We  note  that  the  coefficient  in  this  statement  of  the 
problem  becomes  nonintegrable  (see  also  [1],  [5]). 
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Equations  (1.5)  and  (2.6)  with  corresponding  boundary  conditions 
were  solved  simultaneously  on  an  electronic  computer  for  a  ■  a  in  the 
range  0  £  7  £  8  with  a  maximum  error  of  A  £  l.lcr4  for  values  of  the 
magnitudes  9"  (0),  go  (0)  and  A  <  10"s  for  values  of  magnitude  z*  (0) . 

During  the  calculations  it  was  revealed  that  the  solution  of  a 
first-approximation  equation  is  very  "sensitive"  to  the  solution  of 
zero-approximation  equations.  Tables  1-4  and  Figures  1-5  show  soma 
rounded-off  (to  save  space)  results  of  the  calculations.  The  data  in 
Table  4  was  computed  for  a  =  a  *  1.  The  data  in  Table  5  differs 
somewhat  from  the  corresponding  data  of  Li  and  Nagamatsu  [2],  [5] 
mainly  due  to  the  magnitude  IQ.  For  instance  in  the  latter  study  [3] 
with  H  =  1.4,  values  are  given  of  IQ  =  1.322,  0.9703,  and  0.5956, 
respectively  for  Tw/T0  =  1,  0.6,  and  0.2.  The  assumption  of  the 
presence  of  such  errors  [2],  [3]  was  stated  before  ([7],  p.  359).  In 
Figs.  1-3  the  solid  lines  correspond  to  data  when  T^/To  **  1,  the  dashed 
lines  when  Tw/P0  =  0.6,  the  dot-dash  lines  when  Tw/T0  *  0.2,  £1  »  £». 
The  magnitude  cuj  is  introduced  by  formula 

T/T a  —  0  —  =  g  —  <p'*  -r  0(gi  —  2<p'C')  =  <a0  +  *“»• 

In  the  case  of  an  insulated  plate,  the  slip  increases  the  gas 
velocity  (especially  in  the  near-wall  portion  of  the  boundary  layer) 
and  decreases  the  gas  temper ature  which  is  minimum  in  the  middle  of 
the  boundary  layer  (in  thi3  case  a>t  -  2q>»C*).  A  reduction  of  Tw/P0 
increases  the  role  of  the  near-wall  temperature  Jump.  In  the  near¬ 
wall  part  of  the  boundary  layer  £»  somewhat  decreases,  o>i  sharply 
increases  in  the  main  part  of  the  boundary- layer,  which  leads  to  an 
increase  of  gx .  The  presence  of  extremes  of  these  magnitudes  is 
mainly  because  9*  increases  from  0  to  1  when  X  changes  from  0  to  * 
and  £»  decreases  correspondingly  from  value  of  £  (0)  =  9"  (0)to  a 
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value  (»)  ■  0. 

Slip  and  near-wall  terr$>erature  Jump  decrease  the  pressure ,  dis¬ 
placement  thickness,  local  coefficient  of  friction  drag  and,  espe¬ 
cially,  the  heat-transfer  coefficient  from  gas  to  plate.  This  effect 
decreases  with  a  decrease  of  the  ratio  T^/To* 

Analogous  conclusions  in  the  case  of  an  Insulated  plate  have  been 
obtained  by  means  of  an  integral  form  of  equation  of  momenta  in 
Galkin* s  study  [5]  (velocity  profile  in  boundary  layer  is  linear)  and 
in  Takano's  work  [6]  (velocity  prof  lie  Is  a  f  ourth-degree  polynomial), 
while  In  the  latter  only  the  correction  for  pressure  due  to  slip  when 
x  =  1.4  was  calculated.  In  case  Tw/T0  =  1  when  x  =  1.4,  according  to 
our  data  and  others  [5],  [6]  we  have  respectively:  Pi  *  0.512,  0.20, 
0.515,  &i  =  0.25,  0.15,  =  0.54,  and  0,21.  Hence,  representing 

the  velocity  profile  as  a  fourth-degree  polynomial  Is  sufficiently 
accurate .  A  linear  velocity  profile  can  be  used  only  for  very  simple 
approximate  estimations  of  the  slip  effect,  especially  in  a  case  of 
such  flows  about  insulated  thin  bodies  of  finite  thickness  (for 
instance,  a  wedge),  when  in  the  zero  approximation  motion  is  not  self¬ 
similar. 

The  authors  are  grateful  to  A.  A.  Nikol*skiy,  V.  P.  Shidlovskiy, 
and  M.  N.  Kogan  for  a  discussion  of  the  results  of  this  study. 
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TABLE  1 


x  =»  1,4 


IVT.-I  j 

j  T(tf/T0  -  0.6 

Tw/T9  - 

6.7 

-  i 

1 

) 

9' 

I  U 

—  j 

o.l 

0,075 

1 

1  0,034  j 

0.007 

0.020 

9, 050 

0.050 

0,239 

0.039 

0.2 

0, 147 

O.UOG 

!  0,132 

o,i>;o 

0, 076 

0,116 

0,279 

0.07K 

0,3 

0,210 

i  0,097  j 

j  0,1 95 

O.COl 

0,112 

0,173 

0,318 

0,115 

0,4 

0,282 

!  o,127 

!  0,250 

o.tei 

0, 1  47  | 

0,229 

0,358 

0,  150 

0.5 

0,345 

0,154 

I'.OIC 

0,701 

0, 17  J 

0,28  4 

0,307 

0,184 

0,6 

0.405 

0,178  | 

10,373 

03  721 

0.210 1 

0,3.13 

0,430 

0,210 

0,7 

0,402 

0,201  ! 

'0.427 

0,740 

o.r.s 

0.3ft) 

0,474 

0.245 

0,8 

0,510 

0,220  | 

0, 480  j 

0  700 

0,262 

0,441 

0,512 

0,271 

0,0 

0, 500 

Of 237  i 

0,529 

0,770 

0,281 

0,480 

0,519 

0, 294 

1.0 

0,014 

0, 250 

0,576 

0,797 

0,3' ‘2 

0,536 

0.585 

0.313 

1.2 

0, 008 

0,200 

0,  002 

0,832 

0,326 

0,623 

0,655 

0,340 

1.4 

0,770 

0,208 

0.757 

0,8*15 

0,333 

0,700 

0.719 

0.352 

l.C 

0,829 

o,r>o 

0,600 

0,b.*l 

0,325 

0,703 

0,777 

0.316 

1,8 

0,870 

0,234 

0,652 

0,919 

0,202 

0,824 

0,828 

0,327 

2.0 

0,012 

0, 204 

0,893 

0.950 

0,269 

0,871 

0.871 

0,296 

2.4 

0,900 

0, 136 

0,919 

0.970 

0,187 

0,936 

0,933 

0.216 

2.8 

0,984 

0,076 

0,979 

0.9S7 

0,100 

0,972 

0,970 

0.130 

3.2 

0.994 

0,036 

0,992 

0,995 

o,ar>4 

0,989 

0,988 

0,067 

4,0 

0,999 

0,004 

0,999 

1,000 

0,003 

0,999 

0,999 

0,010 

5,0 

1,000 

0,000 

1,0001 

1,000 

0,000 

1,000 

1,000 

0.000 

TABLE  2 


*=»/. 


A 

rJT‘ 

m\ 

TJt,  =  0.6  | 

rjr.- 

-1 

■ 

9' 

— i 

9' 

-*  s 

r  1 

~ z 

*•  1 

ft 

— J 

0.1 

0,683 

0,033 

0,073 

0,621 

E 

0,040;!  1 ,2 

0,735 

0,231 

0,092 

0,836 

0,2 

0,163 

0,005 

0,144 

o,a;i 

0,078  „  1,4 

0,803 

0,226 

0,704 

0,608 

0.3 

0.238 

0,094 

0,212 

0.C62 

0,114  1,6 

0,857 

0,211 

0,824 

0,897 

0,294 

0.4 

0,120 

0,277 

uxm 

0,148  1,8 

0,599 

0,158 

0,872 

0.2C8 

mam 

0.377 

0,145 

EE  fl 

0,71*3 

0,179  2,0 

0.930 

0,1G0 

0,909 

HRn 

0,234 

0.6 

0,723 

0,208 !  2,4 

0,101 

0  953 

np/fi 

0,157 

KlYB 

0,499 

0,185 

0 ,455 

0,713 

0,233  r  2,8 

0,989 

0,054 

0,9S4 

0,088 

0,088 

0,8 

0,554 

0,201 

0.509 

0,763 

0,255  3,2 

0,99G 

0,024 

0,994 

0.095 

0.042 

0,9 

0,0*5 

0,213 

BEE 

0,782 

0,274  ‘  4,0 

1.000 

0,003 

0,908 

0,005 

1.0 

0,653 

0,222 

^0,607 

0,801 

0,289 ,*  5,0 
fi 

0,000 

1,000 

1,000 

0,000 

TABLE  3 


TJT' 

f'10) 

0 

B 

B 

mm 

B 

St, 

X=:l,4 

1 

0,6 

0,2 

0,01 

0,7627 

0.C775 

0.5S33 

0,5443 

1 

1 ,3003 
0,0407 
0,5502 
0,3728 

0,5091 

0,3657 

0,2174 

0,1449 

B 

vm 

■TrffrM 

0,1531 

0,1150 

0,0926 

V. 

1 

0,6 

0,8544 
0,7429  I 

0 

0,2063 

1,2104 

0.87G4 

0,9089 

0,6532 

0,8527 

0,7229 

0,8146 

0,6004 

0,2084 
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